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Errata 


T:ie  author  has  changed  positions  and  is  now: 


JohnH.  Nath 
R  search  Engineer 
Department  of  Oceanography 
Oregon  State  University 
Corvallis,  Oregon  97331 


Page  19,  Eq.  (26),  should  read: 


Page  20,  the  end  of  the  second  line  below  Eq.  (30),  should  read; 

B  =  2.  /  sin  ^ 

0  a 

Page  the  third  line  from  the  bottom  should  read: 

'•is  a  surface  following  disc  and  the  other  is  a  non  surface 
following,  or" 

Page  33,  the  third  word  In  the  first  line  of  the  last  paragraph  should  be 

"output" 

Figure  42,  in  the  caption  the  depth  should  be: 


"Depth  =  1000  FT.  " 
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ABSTRACT 


Analytical  and  niuierical  studies  were  made  of  single  point 
moorings  of  a  large  disc  buoy  in  deep  water.  Waves  of  different 
frequencies  and  winds  of  different  magnitudes  were  imposed  on  moorings 
of  different  scopes  for  nylon  ropes.  The  numerical  model  results 
compared  favorably  with  analytical  solutions  for  a  straight  vibrating 
string  and  the  buoy  motion  was  validated  with  the  results  from  a 
hydraulic  model  study.  Approximate  calculations,  based  on  the  straight 
vibrating  string  solittion,  were  compared  to  the  numerical  results  of 
curved  mooring  lines  of  various  scopes.  Tne  analytical  prediction  of 
line  tension  was  from  fair  to  poor  but  the  prediction  of  the  position 
of  the  nodes  of  tension  and  the  natural  frequencies  of  the  modes  was 
good.  A  new  method  of  presenting  mooring  line  tension  data  in  dimen¬ 
sionless  fora  is  presented  that  reduces  moorings  of  different  frequen¬ 
cies  and  line  diaaseters  to  a  comanon  scale.  The  versatility  of  the 
numerical  program  was  illustrated  with  two  runs  for  laoorings  in  shallow 
water  with  large  waves  in  addition  to  a  computation  of  line  tension, 
position  and  velocity  for  the  condition  ot  a  wind  shift  of  180  degrees 
on  the  buoy.  The  usefulness  of  the  numerical  piogram  as  a  design 


tool  was  established. 
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INTRODUCTION 


Background  Information 

This  report  will  present  the  current  status  of  research  into  the 
numerical  analysis  of  the  motion  of  the  mooring  line  and  buoy  of  a 
single  point  mooring  in  the  deep  ocean  for  an  oceanographic  buoy.  The 
work  originally  started  in  the  summer  of  1967  for  the  Convair  Division 
of  General  Dynamics  and  was  continued  in  the  summer  of  1968.  During 
the  academic  year  1968" 1969  the  work  was  continued  by  the  author  on  a 
part  time  basis  at  Oregon  State  University  under  the  auspies  of  the 
Office  of  Naval  Research.  Such  a  part  time  arrangement  was  continued 
during  the  past  academic  year,  also  for  the  Office  of  Naval  Research. 
The  objective  of  this  study  has  been  to  develop  a  flexible  numerical 
program  that  can  be  used  for  design  work  and  for  engineering  research 
on  the  aooring  of  an  oceanographic  buoy  in  deep  water  with  a  single 
line.  It  was  desired  that  the  prograR  be  able  to  handle  situations 
with  relatively  large  scopes  as  well  as  taut  iaoorings.  The  effects 
from  waves,  or  problems  on  transients,  such  as  the  anchor  last 
deple>Tscj5t  pro'edure,  can  be  studied.  An  aitcmative,  or  supplessental , 
approach  would  have  been  to  consider  a  solution  ba.sed  on  small 
perturbations  about  the  equilibrium  position  of  the  line.  However, 
such  a  procedure  considers  only  sinusoidal  input  function^  ard  ti-e 
program  that  was  developed  was  such  sore  flexible  in  being  able  to 
cenpute  results  for  a  wide  variety  of  problems  which  ir.clsK’ad  ccttain 
non-linear  effects.  For  design  purposes  it  would  be  it-eful  to  have 
both  t)^«s  of  solutions. 
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The  numerical  program  can  be  considered  to  exist  in  two  parts  - 
one  that  determines  the  static  or  equilibrium  position  and  tension  in 
the  line  due  to  the  action  of  wind,  current  and  gravity  and  the  second 
part  that  determines  the  dynamic  tension  and  motion  of  the  line  due  to 
waves  or  other  forcing  functions.  The  first  part  generates  the  initial 
conditions  for  the  second  part.  Ail  loads  and  system  responses  are 
considered  to  be  acting  in  one  plane. 

The  static  part  of  the  work  has  been  described  in  Refs.  10  and 
11  and  the  dynamic  portion  of  the  program  has  been  described  in  Refs. 

11  and  12. 

The  work  has  been  successful  to  the  point  of  developing  and 
dVouj'ging  a  program  which  predicts  the  position,  velocity  and  tension 
in  ail  parts  of  the  line  and  the  buoy  motion.  The  forcing  function 
can  be  a  wave  of  any  fens,  providing  the  water  particle  velocities 
and  accelerations  are  known  as  functions  of  water  surface  elevation. 

In  addition,  other  transient  problems  can  be  studied,  such  as  an 
anchor  last  deployment  in  still  water,  or  a  transient  horizontal 
motion  of  the  buoy  due  to  a  directional  change  in  the  wind.  In  many 
cases  the  internal  or  structural  dashing  in  the  line  has  been  con¬ 
sidered  as  well  3S  the  hydrodynamic  damping. 

The  differential  equations  of  motion  of  the  line  have  been 
derived  in  Ref.  12,  The  method  of  characteristics  was  used  to  set 
up  the  computational  scheme  for  the  numerical  solution.  Schram  and 
Beyle  (15)  stated  it  well  that  the  method  %:  characteristics  does  not 
liisit  one  to  the  consideration  of  small  perturbations  about  the 
equilibrium.  It  was  also  shown  in  Ref.  12  how  transfer  functions  can 
be  developed  and  a  general  agreement  was  illustrated  between  a 
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computed  transfer  function  and  one  developed  from  experimental  work 
on  a  mooring  in  13,000  feet  of  water  near  Bermuda. 

Scope  of  This  Study 

This  study  continues  tc  consider  the  co-planar  or  two-dimensional 
problem  only.  The  buoy  considered  was  the  large  discus  buoy  (forty 
feet  in  diameter)  developed  by  General  Dynamics  for  the  Office  of 
Naval  Research.  The  mooring  line  characteristics  are  typical  for 
man-made  fibers,  being  non-linear  in  the  stress-strain  relationship. 

An  approximation  to  line  dairoing  was  included.  Some  time  was  required 
for  additional  debugging  of  the  program  from  the  1969  results.  It  is 
felt  that  the  program  is  highly  developed  now  and  can  accomodate  most, 
if  not  all,  configurations  of  interest. 

Favorable  comparisons  were  made  between  the  closed  mathematical 
solutions  for  a  relatively  straight,  vibrating  string  and  the 
numerical  model.  Both  longitudinal  and  transverse  vibrations  were 
considered. 

Some  recent  laboratory  experimentation  with  models  by  General 
Dynamics  has  resulted  in  data  that  enabled  the  author  to  validate  the 
portion  of  the  numerical  model  thac  predicts  the  prototype  buoy 
motion  due  to  different  excitations.  In  addition,  recent  wind  tunnel 
studies  at  Colorado  State  University  have  enabled  a  more  accurate 
assessBient  of  the  wind  di'ag  forces  on  the  buoy. 

Of  some  interest  of  late  has  been  tjie  question  of  how  the 
mooring  line  tension,  and  particularly  that  at  the  anchor,  changes 
when  a  sudden  ISO  degree  wind  shift  occurs  at  the  buoy,  holding  all 
other  variables  constant.  Ihis  problem  was  successfully  studied  and 
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it  was  determined  that  no  unusual  forces  were  exerted  at  the  anchor 
that  were  not  also  distributed  down  the  mooring  line. 

The  primary  interest  in  this  research  was  to  obtain  some  deter¬ 
mination  as  to  the  influence  on  dynamic  line  teixsion  due  to  line 
scope  and  the  magnitude  of  the  wind  as  well  as  the  frequency  of  the 
•vaves.  In  addition,  some  minor  runs  were  made  at  shallow  depths  to 
see  how  the  mooring  line  behaved  when  a  fairly  large  wave  was 
subjected  to  the  system. 

It  is  felt  that  much  more  needs  to  be  done  in  order  to  gain  a 
complete  view  of  the  entire  range  of  problems  associated  with  the 
type  of  moorings  that  were  studied.  However,  due  to  a  real  limitation 
in  time  and  other  resources  it  was  necessary  to  leave  the  very 
interesting  additional  work  to  future  studies. 

Review  of  Recent  Literature 

Reference  12  presents  a  con5)lete  derivation  of  the  equations  of 
motion  of  the  mooring  line  and  buoy.  T5ie  numerical  solvition  of  the 
eqiiations  indicated  that  the  relationship  between  wave  height  and 
mooring  line  tension  was  somewhat  linear  for  a  wave  period  of  7.15 
seconds.  Hence,  transfer  functions,  wh'ch  relate  the  wave  spectrum 
to  the  tension  spectrum  at  a  position  on  the  line  were  developed  for 
a  simulated  mooring  lepiesenting  real  prototype  conditions  near 
Bermuda  in  13,000  feet  of  water.  General  ajrretaent  was  obtained 
between  numerical  predictions  and  field  measured  values  of  the  trans¬ 
fer  function.  In  addition,  frequency  response  curves  for  the  large 
disc  buoy  in  the  heave  and  pitch  modes  of  motion  were  generated  and 
it  was  estiaa.ted  that  the  buoy  was  a  nearly  perfect  surface  follower 
for  all  wave  frequencies  up  to  0.28  cps,  where  the  response  dropped 
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off  sharply.  Waves  of  higher  frequencies  prtjvided  very  little 
stimulation  to  the  buoy  motion. 

It  is  interesting  to  note  that  Devereux,  et.al.,  (2)  reported  that 
line  tension  at  the  buoy  was  quite  linearly  related  to  wave  h  ght  for 
waves  of  mixed  frequencies  for  measurements  taken  at  a  prototype  moor¬ 
ing  in  the  Gulf  Stream  where  the  water  depth  was  about  1040  feet  and 
the  length  of  nylon  line  plus  chain  was  1.330  feet.  The  measurements 
were  taken  during  Hurricane  Betsy  when  the  wave  heights  were  as  large 
as  SO  to  so  feet  and  the  predominant  wave  period  was  ten  seconds. 

The  problem  of  predicting  the  static  position  of  a  buoy  and 
mooring  due  to  steady  current  has  received  considerable  interest. 

Nath  and  Felix  (10)  presented  basic  equations  and  developed  a  numeri¬ 
cal  solution  based  on  equal  increisents  of  tangency  angle.  Bertejujx 
(1)  also  presents  the  basic  equations  and  typical  current  profiles  to 
be  used  off  the  East  Coast  of  the  United  States  and  he  presents  some 
comparative  results  of  calculations  which  study  the  influence  of  line 
diameter  and  scope  on  mooring  design.  .Martin  (7)  describes  ■»  macerical 
program  based  on  equal  line  segments  and  gives  information  on  experi¬ 
mental  drag  forces  on  some  buoy  shapes.  The  report  is  mainly  concerned 
with  mooring  lines  that  have  steel  rope  for  the  upper  portion  to 
resist  fish  bi'e  and  nylon  rope  for  the  lower  portion. 

Treatments  of  the  problem  of  solving  for  the  dynamic  tensions  in 
different  types  of  moorings  have  been  developing  in  recent  years. 

Nath  (12)  describes  the  foundations  for  the  present  study.  Fofonoff 
and  Garrett:  (4)  developed  approximate  theoretical  computations  based 
on  drag  forces  acting  on  the  upper  portion  only  of  taut  moorings. 

They  were  also  concerned  with  the  transient  forces  due  to  the  anchor 
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last  launching  technique.  Millard  (9)  presents  static  tension 
variations  and  particularly  dynamic  tension  variations  wherein  the 
dynanic  tension  variations  were  very  roughly  proportional  to  the  wind 
speed.  He  attributed  that  behavior  to  the  increased  sea  state  for 
higher  winds  but  this  report  will  show  that  it  is  also  reasonable  to 
attribute  at  least  a  large  portion  of  the  increased  dynamic  tension 
to  the  increased  mean  line  tension. 

Goeller  and  Laura  (6)  present  equations  of  motion  for  closed 
solutions  in  complex  fom  for  straight  ccapoimd  lines  where  the  central 
interest  is  in  salvage,  or  raising  loads  from  the  ocean  floor.  They 
included  concepts  of  internal  daisping  in  the  line,  comparing  the 
results  cf  a  sinple  distributed  mass  having  the  usual  two  parameter 
model  of  damping  with  a  lusped  mass  three  parameter  model.  Experiment 
and  theory  showed  that  the  siople  treatment  yielded  good  results  in 
the  region  of  resonance;  where  the  effects  of  dajping  are  most 
iig)ortant. 

Sdtr^  and  Reyle  (J5)  used  the  method  of  characteristics 
transformation  to  the  partial  differential  equations  for  a  flexible 
line  to  develop  a  three-dimensional  numerical  solution  for  the 
dynamic  response  of  cable  towed  systems  at  sea.  They  considered  the 
cable  to  be  incxtensible  and  they  did  not  consider  hydrodynamic  added 
mass.  Since  the  modulus  of  elasticity  of  the  material  was  not 
considered,  the  slope  of  all  characteristic  curves  were  equal  to  the 
velocity  of  a  displacement  wave  on  the  line.  Nath  (12)  showed  that 
the  slope  cf  two  of  the  four  characteristic  curves  were  equal  to  the 
celerity  of  a  longitudinal  elastic  wave  in  the  line.  The  cable 
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towed  systems  illustrated  that  the  transverse  motion  was  damped  to  a 
much  greater  degree  than  the  longitudinal  motion, 

Ihe  anchor  last  deployment  procedure  for  buoy  moorings  was 
investigated  by  Froidevamc  and  Schol:en  (5).  The  analysis  for 
numerical  solution  utilized  a  lumped  mass  representation  of  the  line. 
The  analysis  was  done  for  a  short  system  and  the  results  were  extrap¬ 
olated  for  a  largei  one.  When  the  elasticity  of  the  line  was  intro¬ 
duced  into  the  solution  a  prohibitive  amount  of  computer  tj.ase  was 
encountered, 

A  major  unknown  in  the  analysis  of  buoy  mooring  systems  has  been 
the  hydrodynamic  drag  and  added  mass  characteristics  of  the  buoys, 
Mercier  (8)  describes  an  e,xperimental  model  study  on  various  buoy 
shapes  in  a  wave  basin  and  Felix  (3)  describes  a  model  study  for  the 
large  disc  buoy. 


o 

THEORETICAL  CONSIDERATIONS 

This  section  will  present  a  summary  of  the  final  developments  of 
the  scheme  for  the  numerical  modeling  for  a  single  point  mooring  of 
the  large  ONR  disc  buoy.  The  assumed  line  is  a  typical  nylon  plaited 
rope  with  a  non-linear  stress-strain  diagram  and  simple  two  parameter 
damping,  which  is  inversely  proportional  to  the  frequency,  will  be 
used. 

Governing  Equations  for  Line  Motion 

The  derivations  for  the  governing  equations  for  the  line  motion 
are  given  in  Ref.  (12).  The  method  of  characteristics  was  utilized 
to  transform  the  partial  differential  equations  to  total  differential 
equations.  The  four  dependent  variables  involved  are  the  velocity  of 
the  line  taken  in  the  axial  direction,  which  is  tangent  to  the 
curvature  of  the  line,  v  ,  the  velocity  of  the  line  taken  in  the 

A 

radial  direction,  which  is  perpendicular  to  the  line  tangent,  , 
the  angle  the  line  tangent  makes  with  the  horizontal,  6  ,  and  the 

line  tension,  T  .  The  two  independent  variables  are  the  distance, 
s  ,  measured  along  the  line  from  the  coordinate  axes  as  shown  in 
Fig.  1,  and  the  time,  t  .  The  solution  proceeds  numerically  on  an 
imaginary  s-t  grid  composed  of  eqi;al  time  increments,  At  ,  and 
distance  increments.  As  ,  (except  at  the  upper  boundary  condition, 
which  is  described  later).  Illustrations  and  a  description  of  the 
procedures  are  given  in  Ref.  12. 

The  finite  difference  fora  of  the  total  differential  equations 
can  be  expressed  with  the  following  matrix  equation. 
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On  the 

characteristic 

curve 


Wherein  the  f  teras  represent  the  forcing  functions  on  the 
line.  The  superscript  2  stajids  for  the  values  at  t+&t  and  the 
superscript  1  stands  for  the  values  at  tiae  t  .  That  is,  the 
coefficient  matrix  contains  the  values  of  v  ,  etc.,  at  time  t  . 
In  Eq.  i. 


where  A  is  the  cross-sectional  area  of  the  line,  E  is  the  modulus 
0*^  elasticity,  or  the  slope  of  the  tangent  to  the  stress-strain 
diagram  of  the  line,  y  is  the  saturated  s^ass  density  of  the  line  in 
slugs  per  foot  and  H  is  the  virtual  mass  of  the  line  per  foot  (the 


10 


actual  reass  plus  the  hydrodv'na^a c  added  mass).  For  conpisteness,  the 
forcing  functions  will  be  Jiacsd  below. 


f  =  -  dt  g  (1  ■  ^)  iine  +  V  -  vj  -  “  -  At  (4) 


t*  »  -  At  g  (1  -  sine  +  V  -  V  8  +  v 
2  Oar  C^,.  E 


C  Q  .2 

a  S'^e 


(S) 


fj  »  At 


2 

'r<D 


2  ]  V  |v|  ♦  ( ^ 

Cj  -  I  i  i-  c^, 


\ 


'•  G  -  1 


(6) 


f ,  =  At  r"^^!  1  '^i 

4  [_«D  ■•■  Gi  I  } 


^ "  *^1  w 

p— — j=— I  A  cose  -  A  sine 

^  ^1/  M  ^  ; 


I  Q  .  I  i 

(s-rT^j  s 


V  ♦  (v  +  c  )  e 
r  a  r 


(7) 


In  the  above  equations,  g  is  the  acceleration  of  gravity,  G  ts 
the  specific  gravity  of  the  saturated  line,  Q  is  the  daaping 
coefficient  which  will  be  described  later,  t  is  the  strain  in  the 
line,  Cjj  is  the  drag  coefficient  perpendicular  to  the  line,  D  is 
the  line  diareeter,  Cj  is  the  added  reass  coefficient  (which  is  equal 
to  1.0  for  a  snooth  circular  cylinder),  V  is  the  relative  velocity 
between  the  line  and  the  water  in  the  radial  direction  and  A  and 
are  the  accelerations  of  the  water  particles  in  the  z  and  x 
directions.  The  values  of  the  variables  in  Eqs.  4  through  7  are 
determined  at  time  t  by  an  interpolation  procedure  from  the  s-t 
grid  intersection  points  as  described  in  Ref.  12.  The  second 
derivative  of  the  Im*-  strain  was  determined  simply  with  the  finite 
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Governing  Equations  for  Buoy  Motion 

The  motion  oi,'  the  buoy  was  simply  treated  in  terms  of  rigid  body 
n»tion.  All  the  forces  acting  on  the  buoy  were  considered  in  the 
x2  and  z2  directions  which  ere  shown  in  Fig.  I.  The  acceleration 
of  the  buoy  was  dett rained  at  each  time  station  with  Newton's  second 
law  of  motion  and  the  displacement  and  velocity  for  the  next  tiae 
station  was  determined  by  means  of  recurrence  formulae.  The  forces 
on  the  buoy  were  those  due  to  pressure,  or  buoyancy,  wind  and  current 
drag,  added  mass,  line  tension  and  gravity.  The  forces  were  susased 
in  the  x2  and  z2  directions  and  moments  were  susaued  about  the 
center  of  gravity,  "'‘o  approximate  the  true  di.stribution  of  forces  on 
the  buoy  the  buoy  was  divided  into  a  nuriser  of  pie-shaped  pieces. 

It  was  assumed  chat  dra;'  and  added  mass  forces  were  concentrated 
near  the  lower  chine  of  the  buoy.  By  taking  summation  of  force.«  in 
the  x2  direction  equal  to  the  mass  of  the  buoy  times  the  acceler¬ 
ation  of  the  center  of  gravity  in  the  x2  direction  and  by  subsequent 
re-arrsngment  the  following  equation  can  be  derived. 

KjX2  .  =  Nt  sin>  >  ’  ^x2  *  ^Px2  ^  ^2 

*  ^IR®^x2 


And  in  the  z2  direction. 
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IC3Z2  -  «  cos  ,  .  -  V  •  ‘  Pz2 


(10) 


Likewise,  by  taking  suasnation  of  rooaents  about  the  center  of  gravity 
equal  to  the  laonient  of  inertia  about  the  center  of  gravity  times  the 
angular  accelation. 


KjX2  -  K^z2  ♦  K.^4  =  ^2  *  ^x2  *  ^x2  ^  Ep(soia.  arias) 

■5]''z2  *I  *  ‘1  *X2  -  ‘^1*=”^  ^*22  *  1«  ■  (“) 

< 

In  the  above  equiitior.s  ^  represents  the  derivative  of  ^  with 
respect  to  time  and  4-  is  the  second  derivative.  The  coefficients 
on  the  left  of  the  eqiiaticns  are. 


’i 


CjaP  lY 


^4  '  ’-lA^  ^''^1 


h  "  h 


S  “  '  "  22  Ey  .  CjaP  J:yX2 


(12) 


(13) 


(14) 


(15) 


(16) 

(17) 


The  definitions  of  the  various  terms  in  Eqs.  9  through  17  are: 
ht  is  the  total  buoy  weight,  Is  the  wind  force  cojg)onent  in  the 

x2  direction,  likewise  for  h'^2  »  ^’^2  ^-2  force 

cofifponents ,  p  is  a  pressure  force  on  a  pie-shaped  piece  of  the  buoy, 
the  details  of  which  will  not  be  presented  here,  1^  the  total 

wave  and  current  drag  force  acting  on  the  buoy  in  the  x2  direction 
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and  it  wa5  assumed  that  it  acted  at  a  point  mid-way  between  the 
attachment  point  and  where  the  water  surface  profile  intersects  the 
zl  axis,  which  was  designated  the  distance  Zj  from  the  center  of 
gravity,  is  the  wave  and  current  drag  fores  acting  on  a  pie- 

piece  in  the  s2  direction  based  on  the  relative  velocity  at  the 
IC'-er  chine  and  acting  at  a  distance  from  the  center  of  gravity, 

Daraliel  to  the  x2  axis,  C.-  is  the  inertia  or  added  aass 
coefficient  for  the  buoy  in  the  radial,  or  x2  direction,  t  ,  is 
the  inertia  coefficient  in  the  axial  or  %Z  direction,  o  is  the 
mass  density  of  sea  water,  A_«  is  the  coaponent  of  the  acceleratiors 
of  the  water  particles  at  in  the  x2  direction,  A^„  is  the 

acceleration  of  the  water  particles  at  the  chine  of  a  pie-pisce,  y 
is  the  dispiacesent  volume  of  a  pie^piece,  Z^  is  the  distance  from 
the  center  of  gravity  to  where  the  wind  force  is  assumed  to  be 
concentrated,  2^  is  the  distance  frois  the  center  of  gravity  to  the 
attachment  point,  (swr^.arms)  refers  to  the  various  mosent  arms  from 
the  center  of  gravity  to  the  pressure  forces  on  the  pie-piece  md  q 
is  a  daasjing  coefficient  that  was  used  in  conj‘_.!Ction  vf,th  |  to 
provide  the  proper  damping  and  period  of  the  buoy  as  predicted  by  tne 
nuraericai  model  in  order  that  the  pitch  mode  natch  the  results  of  a 
hydraulic  model  study. 

Boundary  Conditions 

At  the  anchor  end  of  the  line  the  velocities  and  displacements 
are  set  equal  to  tcro  for  ail  time.  Then  the  set  of  four  equations 
and  four  unknowns  as  expressed  in  Eq.  j§  reduce  1:o  two  equations 
and  two  unkrit-wns  for  deterring  the  values  of  9  and  T  at  the 


anchor. 
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The  line  is  divided  into  a  nunsbor  of  equal  segments,  As  long, 
at  time  =  0.  The  segment  lengths  reusain  fixed  in  time  except  for  the 
end  segment  which  attaches  to  the  buoy,  the  length  of  which  is 
determined  at  each  time  increment.  Tlte  nusfoer  of  segments  can  be 
increased  or  decreased  according  to  whether  the  line  is  lengthening 
or  shortening.  The  procedure  is  to  use  the  method  of  characteristics 
to  determine  v  ,  v  ,  b  and  T  for  all  segments  except  the 
one  next  to  the  buoy.  Knowing  As  and  i  for  each  segment,  the 
coordinates  at  the  ends  of  each  segment  are  determined.  Since  the 
coordinates  of  the  attachment  point  at  the  b-oy  are  known  from  the 
buoy  action  subroutine,  the  length  of  the  last  segment  can  be  deter- 
mined,  which  also  determines  6  for  the  last  segment.  Knowing  the 
distribution  of  tension  in  the  line,  tue  strain  in  each  equal  segmen*-. 
can  be  detenained  and  since  the  total  length  of  the  line  is  known  the 
strain  in  the  end  segment  can  be  calculated,  thence  the  end  tension 
from  the  stress-strain  diagram.  The  velocities  of  the  end  segment 
are  determined  from  the  buoy  stibroutine.  ^w^e  details  on  the  above 
procedures  ar*  presented  in  Ref.  12. 
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PHYSICAL  CHAPAa'ERISTICS  OF  THE  SYSTEM 

This  study  considered  rcooring  lii es  which  are  characterized  by 
plaited  nylon.  Deep  water  depths  were  mostly  considered  because 
several  harTRonics  ut  the  resonant  frequencies  for  such  mooring  con¬ 
ditions  can  exist.  Line  diameters  of  2.0,  2.5  and  5.5  inches  were 
used  with  scopes  of  0.80,  0.87  and  1.18,  where  scope  is  the  original 
unstreteched  length  of  the  line  divided  by  the  water  depth.  Two 
conditions  of  wind  and  c  rent  were  investigated.  The  winds  were  50 
knots  and  ISO  knots  and  the  corresponding  current  profiles  that  were 
assumed  are  shown  in  Fig.  2. 

It  is  very  difficult  to  determine  the  elongation  characteristics 
for  ropas  made  of  synthetic  fibers  for  a  general  study.  Generally, 
pciycster  materials,  have  less  elongation  than  nylon  for  the  same 
stress  condition.  However,  the  type  of  rope,  twist- lay,  plaited, 
etc.,  and  the  tightness  of  the  strand  will  influence  the  load- 
elongation  curves  to  a  great  degree,  so  that  it  is  possible  and  it 
has  occurred  that  certain  types  of  nylon  ropes  will  display  less 
elongation  than  certain  types  of  polyester  ropes.  Generally,  of 
course,  the  reverse  is  true.  In  addition,  the  stress-strain  charac¬ 
teristics  are  considerably  influenced  by  subo?rsion  of  the  rope,  yet 
practically  no  information  exists  on  this  topic.  Thus,  a  designer 
has  little  information  to  work  with  urless  tests  can  be  performed, 
under  water,  on  the  very  rope  he  plans  to  use.  For  this  study  is  was 
assumed  that  the  mooring  rope  was  subject  to  an  initial  .and  permanent 
strain  of  O.i  from  the  initial,  or  some  subsequent,  loading.  Tue 
stress-strain  diagram  is  presented  in  Fig.  5.  .Also  on  Fig.  3  ar<*  the 
results  of  some  testing  that  are  reported  on  in  Martin  (")  snd  Wilson 
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(’S)  and  the  stress-strain  diagram  used  in  Ref.  12.  The  final 
sclectio.i  of  the  stress-strain  diagram  was  somewhat  influenced  by  test 
information  provided  by  the  Columbian  Rops  Company. 

The  internal  hysteretic  damping  characteristics  for  the  line  were 
estimated  from  Information  given  in  Refs.  15  and  16.  The  procedure 
was  to  measure  the  area  (thus  the  energy  dissipated  )  of  several 
hysteresis  loops,  approximate  them  as  elipses  and  utilize  the  follow¬ 
ing  equations.  For  an  eliptical  hysteresis  loop,  the  an5)litude  of  the 
loading  function  ,  is  related  to  the  amplitude  of  the  displacement, 

,  and  the  phase  shift,  B  ,  with: 


sin  B 


Area  within  the  el ipse 

v  P  X 
o  o 


(18) 


Thus,  B  can  be  determined.  Simple  linear  dajiqjing  in  a  continuous 
system  is  characterized  by 

a(s,t)  =  E  c(s,t)  Q  1“  (19) 


or  in  a  non-linear  system  (as  considered  in  this  report)  by 
o(s,t)  =  o(e)  ♦  Q  If' 

0  L 


(20) 


where  o  is  the  stress  (calculated  as  the  line  tension  divided  by  the 
original  cross-sectional  area  of  the  line).  E  is  the  modulus  of 
elasticity  of  the  material,  e  is  the  strain  and  Q  is  the  dan5>ing 
coefficient.  The  damping  coefficient,  Q  ,  can  be  estimated  by 
considering  the  linear  system,  wherein 

tan  S  =  ^ 

where  u  is  the  frequency  of  the  forcing  function  and  is  the 


(21) 
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mean,  or  dyn.aaic  modulus  ot  elasticity.  Tlius 


Q  =  —  Ej  tan  sin  ^ 

ji!  d 


Thus  from  Refs.  15  anf-’  i6. 


Q  was  estimated  to  be 


Q  = 


lb  -  sec 
ft2 


(22) 


(23) 


For  The  solution  of  the  numerical  program  the  drag  forces  in  the 
longitudinal  direction  on  the  line  were  igncred.  It  has  been  shown 
in  several  studies  that  they  have  negligible  influence  when  scope  is 
less  than  2.  The  drag  coefficient  in  the  radial  direction  was  taken 
equal  to  1.4  and  the  added  mass  coefficient  was  given  the  theoretical 
value  of  1.0  for  a  smooth  circular  cylinder.  The  saturated  mass 
density  of  the  line  was  developed  as 

V  =  1.71  d2  (24) 


where  u  is  the  mass  per  foot  and  D  is  the  line  diameter  in  feet. 
Eqxiation  24  was  based  on  a  saturated  isass  density  of  the  line  of  2.17 
slugs  per  cubic  foot. 

A  physical  description  of  the  buoy  is  presented  in  Ref.  12. 
Several  coefficients  were  derived  for  the  buoy  from  experimental 
information.  Ihe  drag  coefficient  in  the  x2  direction  i.s  a  function 
of  Froude  niasber,  but  for  the  relative  velocities  considered  here  it 
was  felt  it  could  be  assu3sed  to  be  constant  at  0.035.  The  added  mass 
coefficient,  Cjj^  ,  was  0.4. 

In  order  to  obtain  similarity  between  the  numericcl  sjodel  and  tht*. 


results  of  the  hydraulic  model  study  reported  in  Ref.  3  with  respect 
to  response  decay  curves,  it  was  necessary  to  arbitrarily  nxsdify  the 


J8 


usubI  expres'  ion  for  drag  force  in  the  z2  direction  of  the  buoy  '  y 
making  the  force  proportional  to  the  velocity  instead  of  the  velocity 
squared.  In  Eq.  10,  was  established  for  one  pie-piece  as 


^(40)-  1  P  y 
4  12  2  22 


(25) 


Thus  the  dreg  coefficient  is  no  longer  dimensionless.  For  this  work 
the  value  of  was  12.0.  The  value  when  the  relative  acceleration 

in  the  22  direction  at  the  chine  was  positive  of  was  taken  as 

3.0,  and  0.8  when  that  relative  acceleration  was  negative.  The 
auxiliary  damping  coefficient,  q  ,  in  Eq.  11  was  given  the  value 


200,000  ft-lb-sec.  The  distance  X  ,  was  determined  to  be  9.0.  The 


dimensional  wind  drag  coefficient,  C_  ,  where  W  =  Cr.  p  .  V-/2  , 

®  ’  Dw  '  Dw  air  w  ’ 

was  defsra  ned  experimentally  in  Nath  (13)  to  be  140.  The  resulting 


buoy  nation  due  to  the  above  selection  of  coefficients  is  given  in  the 


next  section. 


NUMERICAL  RESULTS 


Longitudinal  and  Transverse  Motion  of  the  Line 

On**  way  to  check  a  numerical  snodel  such  as  the  one  presented  here 
is  to  compare  the  results  from  it  to  those  from  a  closed  solution  of 
a  simplified  problem.  This  was  done  for  the  case  of  a  taut  vibrating 
string,  the  solutions  of  which  are  well  known.  The  development  of 
the  analytical  solution  will  be  presented  first,  followed  b>  coE5)ar- 
isons  of  analytical  results  with  those  from  the  mjmerical  model. 

Consider  a  straight,  weightless  line  held  vertically  as  shown  in 
Fig.  4.  The  line  has  a  cross-sectional  area  of  A  and  a  modulus  of 
elasticity  of  E  .  The  initial  tension  in  the  line  is  T^  .  At 
z  =  L  the  line  is  forced  with  a  longitudinal  excitation,  sin  ut  , 
where  u  is  the  radian  frequency  of  the  excitation,  and  it  is 
assumed  that  the  initial  ten.sion  is  adequate  to  preserve  the  condition 
of  positive  tension  in  the  line  at  all  locations  and  times.  It  is 
assumed  that  the  effects  of  danroing  are  negligible.  The  vertical 
disp]£.»ei5ent  of  a  particle  of  the  line  will  be  designated  as  c  and 
the  strain  in  the  line-  is  then  ,  wliere  the  subscript  indicates 

the  partial  derivative  with  respect  to  z  .  Only  the  steady  state 
vibration  solution  to  the  problem  is  desired.  It  is  well  known  that 
the  governing  equation  is  the  wave  equation,  which  is: 


(26) 


where  is  the  second  derivative  of  the  displacement  with  respect 


to  time  and  a  lums  out  to  be  the  celerity  of  an  elastic  wave  along 
the  line  also  given  by  Eq.  2. 
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The  displacement  is  then  a  function  of  the  coordinate,  z  ,  and 
the  time,  t  ,  and  the  boundary  conditions  are: 


4(0, t)  =  0 

(27) 

?(L,t)  =  Z^  sin  ojt 

(28) 

It  should  be  noted  that  L  >>  and  that  the  problem  has  been 

linearized  by  iicposing  the  boundary  condition  at  z  =  L  instead  of 

Che  true  condition,  z  =«  L  Z  sin  ut  .  The  solution  is  obtained 

0 

by  assuming  that  it  is  a  product  in  the  form: 


4(z,t)  =  {A  cos  ^  z  +  B  sin  —  z)  sin  ut  .  (30) 

By  impofing  the  first  boundary  condition  it  is  seen  that  A  =  0  . 

By  imposing  the  second  boundary  condition  it  is  seen  that  B  =  Z  /  ^ 

0  a 

Tne  solution  for  the  displacement  due  to  dynamic  motion  is  then 


s  r* 

sin  — 
a 


and  for  the  strain: 


,  ,  _  u  1  wr  . 

C.(Zs-)  =  r - ="  cos  —  sin  ut 

z  0  a  .  ojc  a 

a 


Most  ropes  do  not  have  linear  elastic  characteristics,  but  if 
they  are  stressed  to  at  least  20%  of  their  ultimate  strength  the 
stress  strain  diagram  becomes  approximately  linear.  The  linear 


condition  was  assinaed  for  this  section  of  the  work. 


For  a  linear  elastic  material  the  strain  is  directly  proportional 
to  the  stress  through  the  modulus  of  elasticity,  so  that  the  solution 
for  the  dynamic  line  tension,  T  ,  is 


2i 


T(z,t)  =  AE  Z  - 
o  a 


t:-  cos 


sin 


ax 


(i)Z 

a 


sin  u)t 


(33) 


For  the  total  tension  the  steady  tension,  ,  must  be  added  to 

o 

Eq.  (33).  By  substituting  Eq.  33  into  Eqs.  26  through  28  it  is  seen 
that  it  is  the  solution  desired  for  the  steady  state  vibration 
condition. 

Equation  33  shows  that  the  line  tension  has  infinite  values  where 

sin  —  =  0  .  This  will  occur  when 

a 


a>L 

a 


n  IT 


n  =  1,2,3,  . 


(34) 


(The  case  for  n  =  0  is  a  trivial  case  where  either  L  =  0  or 
0)  =  0) 

The  numerical  model  was  tested  by  modifying  it  to  consider  a 
linear  elastic  material  with  no  daE5)ing  and  with  longitudinal  sswtion 
only,  A  frequency  of  0.7  radians  per  second  was  used  to  represent 
ocean  wave  frequencies  that  can  occur  as  strong  swell.  5®veral  water 
depths,  or  lengths  of  line,  were  checked  and  Table  i  shows  the  results 

The  line  tensions  were  plotted  by  the  computer  for  the  anchor 
and  the  buoy  end  as  functions  of  time.  E:-jaraples  of  suen  plots  are 
shown  as  Figs.  5  and  6.  It  was  also  seen  by  inspection  ■:>f  the 
printed  output  that  nodes  of  zero  tension  fluctuation  occurred  at  the 
positions  on  the  line  which  are  predicted  by  Eq.  33. 

It  was  attempted  to  create  a  resonant  condition  by  selecting  a 
frequency  such  that  wL/A  =  2Tt  ,  The  purpose  was  to  investigate  the 
effects  of  internal  damping  on  preventing  infinite  response  as 
predicted  by  Ec  33.  However,  even  by  setting  the  dancing  coefficient 
Q  ,  equal  to  zero,  infinite  response  was  not  obtained,  which 
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indicated  that  the  numerical  solution  as  presented  was  naturally- 
damped.  At  this  writi.''^'  it  has  not  been  determined  why  the  numerical 
solution  as  presented  here  is  naturally  damped. 

It  was  also  discovered  that  the  full  amount  of  damping,  as 
predicted  by  Eq.  23,  produced  a  numerical  instability  in  the  program. 
However,  when  one-tenth  of  the  theoretical  value  of  damping  was  used, 
the  solution  was  stable  for  roost  conditions  and  it  is  felt  ?hat  this 
degree  of  damping  in  conjunction  with  the  natural  damping  will  produce 
results  that  will  approximate  the  action  of  a  real  nylon  line.  The 
results  of  testing  the  solution  for  a  20,000  feet  long  line  are  given 
in  Table  2. 

The  results  show  that  for  the  conditions  assumed  the  oscillations 
in  tension  were  reduced  by  one-third.  For  frequencies  that  were  not 
near  the  resonant  frequencies  the  dampirsg  had  little  influence  on  the 
oscillations,  as  expected. 

The  numerical  model  was  also  checked  in  the  transverse,  or  x 
direction.  In  this  case  the  steady  state  transverse  vibration  was 
determined  for  all  positions,  z  ,  and  time,  t  ,  The  analysis 
assumed  that  displacements  were  very  small  and  that  the  tension  in 
the  line  was  constant  in  time  and  position. 

The  assumed  conditions  were  that  both  ends  of  the  line  remain 
fixed  and  at  time  equal  to  zero  the  transverse  velocity  of  the  line 
was  sinusoidally  distributed  along  the  line.  As  time  progressed  the 
horizontal  displacemei.r  of  the  line  was  also  longitudinally  distributed 
and  the  following  development  will  show  the  closed  solution  for  the 
horizontal  displacement,  n(z,t)  .  The  boundary  conditions  for  this 
problem  are. 
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=  n(L,t)  »  0 


(35) 


The  initial  condition  is 


\  (2,0)  =  sin  g 


uz 


(36) 


where  is  any  convenient  velocity  amplitude  .vid  b  is  the 
longitudinal  celerity  of  a  displacement  wave  al-. ng  the  line  given 
also  by  cq.  3. 

The  governing  equation  is  again  the  wave  equation, 

'tt  tt 

A  product  solution  is  assumed  and  the  general  form  is 
o(2,t)  =  (C  cos  ^  2  ♦  D  sin  ^  z)  sin  oit 


(37) 


(38) 


After  applying  the  boundary  conditions  and  the  initial  condition,  the 
solution  is  found  to  be. 


n(z,t) 


V 

o  .  “z  . 

—  sin  R—  sin  at 

ii)  ^ 


(39) 


Equation  39  shows  that  there  are  nodes  on  the  line  of  zero  displace¬ 
ment  for  all  time.  In  the  numerical  program  the  tensions  were  set  co 
be  constant  with  time  and  position  and  all  drag  and  added  mass 
coefficients  were  set  equal  to  zero. 

For  the  numerical  work,  an  initial  tension  of  3905.2  pounds  was 
assumed  in  addition  to  a  frequency,  ^  ,  of  0.251  rad/sec  and  a  line 

length  of  20,000  feet.  It  is  seen,  then,  from  Eq.  39  that  nodes 
should  occur  in  the  line  at  z  =  5,000,  10,000,  15, 000, and  20,000 
feet.  The  results  cf  the  numerical  work  showed  that  the  nodes 
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occurred  at  about  z  equal  to  5400,  10,000,  15,400  and,  of  course, 
20,000  feet. 

The  maxirauai  magnitudes  of  the  displacements  should  be  18.5  feet, 
as  predicted  by  Eq.  39.  Figure  7  shows  that  the  displacements  were 
not  perfectly  symmetrical,  mostly  due  to  the  fact  that  the  nodes  did 
not  occur  exactly  where  they  were  supposed  to  occur.  However,  the 
average  Maximum  displacement  as  determined  from  Fig.  7  is  18.3  feet. 

The  conclusion  from  the  above  tests  was  that  the  rating  of  the 
nmnericai  model  can  be  described  as  from  good  to  excellent,  the  first 
rating  applied  to  the  transverse  motion  and  the  second  referring  to 
the  longitudinal  motion. 

.Many  buoy  moorings  are  e-stablishes.  with  a  ver>'  taut  line  running 
almost  vertically  from  the  anchor  to  the  buoy.  It  will  be  shown 
below  that  the  resonant  length  or  frequency  of  the  mooring  line  will 
depend  a  great  deal  on  the  type  of  buoy  that  is  used  at  the  water 
surface. 

One  condition  for  consideration  is  that  shown  by  Fig.  4  where 
the  longitudinal  displacensent  at  the  top  of  the  line  is  imposed  as 
the  boundary  condition.  The  equation  of  motion  in  terms  of  displace¬ 
ment  is  given  by  Sq.  31  and  the  line  tension  is  given  by  Eq.  33. 

These  equations  show  that  resonant  conditions  exist  when  wL/a  =  m 

Consider  ne;;t  equal  coraditions  except  that  the  boundary  condition 
at  the  top  of  the  line  is  a  force  given  by 

The  new  boundary  conditions,  for  a  ssaterial  with  a  linear  stress- 
strain  relationship,  are  given  by  Eq.  40  and 
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C^CL.t)  =  ~  sin  ut  (41) 

On  substitution  of  the  boundar)'  conditions  into  Eq.  30,  one  again 
fiiVJs  that  A  =  0  and  that 

B  —  cos  —  sin  iDt  =  -Tp  F  sin  ut  (42) 

a  a  AE  o  ' 


inus. 


B  =  -£  i  1 
A£  u  (i)L 
cos  — 


and  the  solution  for  displacements  is, 
F 


C(2,t)  =  - p 

•  '  AE  -jj  uL 

cos  - 

a 


CJ2 

sin  —  sin  tijt 
a 


(43) 


(44) 


The  curious  difference  between  Eqs.  44  and  32,  is  that  infinite 
responses  occur  for  Eq.  32  at  tuL/a  =  nrr  but  for  Eq.  44  they  occur  at 
uL/a  =  n  y  . 

For  design  purposes  the  resonant  conditions  should  be  avoided. 

For  very  taut  moorings  that  approximate  the  configuration  shown  in 
Fig.  4  the  most  serious  mode  of  motion  will  probably  be  the  longitudi¬ 
nal  one.  Dandling  may  be  relatively  low.  The  resonant  length  of  line, 
given  a  particular  wave  frequency,  will  then  depend  on  the  bovmdary 
condition  at  the  top  of  the  line. 

Consider  two  types  of  buoys  that  may  conceivably  be  ssrored  to 
nearly  vertical  mooring  lines  at  the  air-water  interface.  One  buoy 
is  a  surface  following  large  disc  and  the  other  is  a  non  surface,  or 
relatively  stable,  spar.  For  a  flexible  mooring  line  the  large  disc 
will  iB5>ose  a  displacement  at  the  upper  end  of  the  line  as  it  closely 
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fo] lows  the  undulations  of  the  water  surface  whereas  the  spar  buoy 

will  impose  a  time  varying  force  on  the  line  as  the  waves  pass  the 

spar.  For  a  particular  wave  frequency  the  resonant  length  of  line  for 
the  disc  mooring  will  be 

L  =  n  ^  n  =  1,  2,  2,  .  .  .  (45) 

For  the  spar  buoy,  the  resonant  length  of  mooring  line  will  be 

L  =  n  ^  n  =  1.  3,  5.  .  .  .  (46) 

For  other  buoy  types,  such  as  the  aid  to  navigation  buoys  used  by  the 
U.  3.  Coast  Guard,  some  condition  between  Eqs.  45  and  46  will  exist. 

In  addition,  several  wave  frequencies  exist  at  sea  and  all  resonant 
ccndAtions  with  respect  to  the  entire  wave  frequency  spectrum  should 
be  investigated. 

To  illustrate  with  numbers  the  foregoing  concepts,  consider  a 
common  ocean  wave  period  of  seven  seconds.  The  frequency  will  be 
0.90  radians  per  second.  Assume  that  a  one-inch  diameter  nylon  line 
(that  is  nearly  neutrally  buoyant)  is  used  and  that  it  is  pre-stressed 
to  2200  pounds.  The  modulus  of  elasticity  will  be  about  7.5  x  10^  psf 
according  to  Fig.  3.  The  saturated  density  of  the  line  will  he  about 
0.0118  slugs  per  foot.  Thus  the  celerity  of  a  longitudinal  elastic 
wave  in  the  line  will  be  about, 


a  = 


.00547  X  7.5  X  10^ 
0.0118 


(47) 


or, 


a  *  1860  fps 


(48) 
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For  the  disc  buoy,  the  resonant  lengths  of  mooring  line  will  be 
6,500  feet,  13,000  feet,  etc.  For  the  spar  buoy,  the  resonant  lengths 
of  line  will  be  3,250  feet,  9,750  feet,  etc.  It  should  be  stated 
again  for  emphasis  that  nil  wave  frequencies  should  be  examined  for 
the  resonant  lengths  of  the  line  and  the  type  of  buoy  must  be  given 
careful  consideration. 

Now  consider  the  equations  for  line  tension  for  both  the  above 
conditions.  For  the  disc  buoy  the  tension  is  given  by  Eq.  33.  For 
the  spar  buoy  the  tension  will  be, 

cos  — 
a 

Fcr  both  conditions  the  dynamic  line  ten-Sion  is  a  maximian  at  the 
anchor!  For  the  disc  buoy  the  dynamic  tension  at  the  buoy  may  be  zero 
if  a)L/a  s  nn/2  with  n  -  1,3,5,  ...  .  For  the  spar  buoy  the  d>Tiamic 
line  tension  at  the  buoy  will  be  sin  ut  and  the  displacement  will 
be  zero  if  uL/a  =  ,  n  =  l,2 . 

Buoy  Motion 

A  necessary  part  of  this  study  was  to  select  the  proper  values 
of  the  coefficients  in  Eqs.  9  through  17,  The  procedure  was  to  rely 
heavily  on  the  1:10  scale  model  study  in  Ref.  3.  In  particular  some 
response  decay  curves  for  the  pitch  and  heave  motions  were  generated 
that  iiqilicitly  display  the  effects  of  all  the  coefficients.  By  trial 
and  error  the  coefficients  were  evaluated  for  the  numerical  model 
until  the  predicted  response  curves  nearly  matched  the  best  estimate 
for  the  response  curves  of  the  prototype.  The  response  curves  for 
the  prototype  were  estimated  fro-i  the  response  curves  from  the  scale 
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model  study  by  considering  the  Froude  modeling  relationships  and  the 
shift  in  natural  periods  due  to  different  degrees  of  damping. 
Hydrodynamic  damping  in  the  model  studies  should  be  greater  than  in 
the  prototype  because  of  the  Rejniolds  number  scale  effect.  Tlie 
natural  period  in  heave  for  the  prototype  predicted  by  the  model  study 
was  about  3.6  seconds.  After  accounting  for  the  Reynolds  number 
effect,  the  period  was  reduced  to  3.2  seconds.  Tne  period  in  heave 
from  the  numerical  model  was  3.1  seconds,  which  was  considered  to  be 
close  enough-  Peaks  in  tension  spectra  from  prototype  measurements 
occur  at  periods  of  from  3.1  to  3.4  seconds  and  it  is  predicted  here 
that  they  were  due  to  the  heaving  motion  of  the  buoy.  The  model  study 
also  showed  a  ratio  of  1.128  between  the  heave  period  and  the  pitch 
period.  The  ratio  for  the  numerical  work  is  1.148. 

Figure  8  presents  the  normalized  pitch  response  decay  curve  for 
the  numerical  work  and  the  results  from  the  hydraulic  :./odel  study 
where  the  rscrmalizing  pit^'h  is  the  pitch  at  time  equal  to  zero  and  the 
normalizing  time  is  the  natural  period.  A  similar  curve  is  presented 
in  Fig.  9  for  the  heave  motion;  however,  the  beginning  portion  of  the 
data  from  the  hydraulic  study  appeared  to  be  influenced  by  the  release 
of  the  buoy  model  at  the  start  of  testing.  Thus  the  curve  for  the 
hydraulic  model  was  arbitrarily  shifted  to  the  left  so  that  the 
subsequent  vibration  would  be  in  phase  with  those  from  t'.e  numerical 
work. 

One  minor  but  interesting  phase  of  this  work  was  to  observe  the 
predicted  motion  of  the  buoy-  on  the  surface  of  a  wave.  For  thio  .  -rt 
of  the  work  it  was  assumed  thit  the  buoy  was  moored  in  200  feet  of 
water  with  an  elastic  "-attrial  that  produced  an  eq-ui valent  spring 
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constant  at  the  buoy  of  10  Ibs/ft,  always  directed  toward  the  anchor, 
with  an  initial  tension  of  £0,?00  lbs.  As  in  all  vhe  wor.<  involving 
waves,  the  wave  heigh  %  built  up  linearly  with  respect  to  time  over 
one  wave  period.  Thu.s  Fig.  10  shows  how  the  wave  height  was  increased 
over  the  first  wave  period  atd  how  the  buoy  and  mooring  behaved  for  a 
wave  260  feet  long  and  36  feet  high.  For  this  case  the  mass-spring 
system  in  conjunction  with  the  water  particle  acceleration  was  such 
that  the  attachment  point  was  eicposed  to  the  •'.ir  shortly  after  14 
seconds  and  the  program  was  written  to  stop  foe  such  an  occurrence. 
Figure  li  shows  the  buoy  successfully  negotiating  the  same  wave  for 
the  free-floating  condition  and  Fig.  12  shows  successful  transit 
of  the  moored  condition  when  the  wavs  height  was  reduced  to  30  feet. 

A  much  longer  and  higher  wave  is  press.Tted  in  Fig.  13.  Figure  13 
clearly  shows  the  buoy  heaving  on  the  •  ave  at  the  nati”  4I  damped 
period  of  3.1  seconds,  as  well  as  at  the  period  cf  the  wave. 

Water  Sheave 

Another  minor  but  very  interesting  investigation  made  during  thi.s 
study  was  that  which  is  sometimes  referred  to  as  tho  'Sater  she  /£•'’ 
effect.  An  illustration  of  the  problem  is.  say,  a  mooring  in  20,000 
feet  of  v/ater  subjected  to  a  150  kt.  wind.  At  some  tis.  the  wind 
suddenly  shifts  180°  and  blows  the  buoy  bach  alorig  the  alignment  of 
the  mooring  line  until  it  and  tne  line  reach  new  equilibrium  positi'.n? 
How  the  line  tension  varies  durxng  the  transient  motion  has  been  "r- 
question.  It  is  felt  that  the  form  of  the  line  tension  rolutior 
presented  here,  based  on  the  method  of  characteristics;  is  particular! 
well  suited  to  h.andle  this  problem  in  a  nxnnericai  rcan^'or. 
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\  siH^lar  problem  as  that  posed  above  was  iavesii gated.  Hov  er, 
in  order  to  save  computer  time  the  top  of  the  line  was  transiated  at  i 
higher  velcrity  than  would  bs  the  case  if  the  wind  drag  only  on  the 
biioy  w  TC  re-'ersed  First  it  was  determined  that  if  the  bucy  broke 
free  froor  mooring,  a  150  kt  wind  would  pusn  it  through  the  water 
with  a  speed  of  13  fps.  It  was  decided  to  translate  the  line  at  the 
ccnservativsl.T  much  higher  velocit)  of  25  fps,  increasing  it  linearly 
fren  h  fps  in  the  period  of  20  seconds.  Ar,  initial  "in  place"'  or 
tensioned  scope  of  about  1.3  was  assumed  and  the  initial  r.on- 
equtlibriura  tension  of  7000  lc..s.  was  established  with  a  line  di3«^ter 
of  1.5  inches  for  a  .cateriai  with  ciaracteri sties  as  shewn  in  Fig.  3. 
Tne  tirst  computation  consisted  of  fixing  the  tep  of  che  line  in  place 
and  allowing  it  to  reach  an  equilibrius  position.  oth  positions  are 
shewn  at  time  equal  to  cere  o'.  fig.  Ihen  the  line  was  moved  to 

the  left  with  the  velocity  function  presented  above,  fhe  sequential 

positions  of  the  line  as  detemiiiied  by  the  numerical  ccaputaticr.s 
are  shown  in  Fig.  14. 

The  tension  in  the  line  dropped  during  the  first  150  seconds  and 
then  gradually  increased  until  nearly  rupturing  the  line,  at  wliich 
time  the  cemtu^ations  were  stopped.  During  the  increase  in  tension 
phase  the  aaximum  tension  occurrefl  at  about  0.067  the  distance  down 
the  line.  The  line  length  changed  in  accordance  with  the  state  of 
''ension  along  the  line.  The  tension  contours  on  the  s-t  plane  are 

shown  in  Fig.  IS.  It  should  be  noted  when  examining  Fig.  15  that 

station  1  represents  the  anchoi  and  subsequent  51 ation  numbers  are 
spaced  S6C.7  feet  apart.  It  can  be  seen  that  the  tension  at  the 
anc!.or  decreased  for  about  the  first  150  seconds  and  t)ie>"  steadily 
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increased.  Toward  the  end  of  the  coasputations  the  tension  in  the  line 
was  nearly  constant  with  respect  to  line  position,  which  •$  reasonable 
by  inspection  of  Fig.  14. 

If  the  top  boundary  condition  hau  been  the  actual  buoy  with  the 
wind  onlj  driving  it,  the  tensions  would  have  been  lower  because  cf  the 
lower  buoy  velocity  and  because  af  the  differences  in  the  initial 
loading  coiidition.  When  the  horizontal  cooponent  of  the  line  tension 
at  the  buoy  had  equaled  the  wind  drag  of  10,800  pounds  the  buoy-motion 
would  stop.  Thus  it  can  be  seen  that  the  conditions  imposed  were  much 
more  severe  than  what  would  have  been  experienced  at  sea.  If  the 
submerged  weight  of  the  anchor  were  less  chan  about  10,000  pounds  it 
wo\:ld  have  been  lifted  from  the  bottom  (assuming  no  adhesion  tc  the 
bottom)  at  around  1000  seconds  or  more.  If  the  real  boundary  condition 
had  been  imposed  at  the  buoy  end,  the  motion  would  have  been  quite 
slow  and  the  computer  time  necessary  to  complete  the  problem  would 
have  been  prohibitive  within  the  budget  for  this  study. 

Effects  From  the  Environment  at  Sea 

The  main  purpose  of  this  reseai'ch  was  two- fold;  1.)  to  complete 
the  debugging  of  the  numerical  program  and  to  test  it  against  closed 
mathematical  solutions,  as  already  discassed,  and  2.)  to  gain  some 
insight  into  the  influence  on  dynamic  line  tensions  from  different 
wind  and  current  loads  and  ux.ferent  line  scopes.  The  work  related 
to  the  first  item  has  been  presented.  For  the  second  purpose  it  was 
decided  to  design  a  few  moorings  in  20,000  feet  of  water  for  certain 
scopes  for  a  maximum  wind  condition  of  ISO  kts.,  subjecting  each 
mooring  to  waves  of  various  frequencies.  Tlien  the  same  moorings  ■'ere 
tested  with  the  same  waves  but  with  a  moderate  wind  velocity  of  50  kts. 
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The  current  velocity  profiles  corresponding  to  the  two  wind  velocities 
are  presented  in  Fig.  2. 

As  presented  previously,  the  computer  program  first  determined 
the  steady  state,  or  equilibrium,  position  of  the  lines.  The  f^'e 
different  conditions  are  shown  in  Fig.  16.  In  all  cases  the  equilib¬ 
rium  tension  in  the  line  was  kept  below  47%  of  the  breaking  strength 
as  predicted  by  Fig.  3.  Thus  three  diameters  were  selected,  as  jwn 
in  Fig.  16.  A  summary  of  the  lead  and  scope  conditions  is  given  in 
Table  3, 

Two  short  runs  were  made  in  shallow  water  in  order  to  obttin  a 
brief  view  of  the  influence  of  fairly  large  v^aves  on  line  tenf.ion. 

The  steady  load  conditions  for  these  two  runs  are  summarized  in  Table 
3.  The  profile  equilibrium  positions  have  not  been  included  in  a 
Figure  because  they  were  practically  straight  lines. 

The  dynamic  program  based  on  the  method  of  characteristics  was 
compared  to  the  steady  state  program  in  terms  of  the  static  loads. 

The  procedure  was  to  assume  a  water  depth  of  10,100  feet,  a  scope  of 
2.02,  diameter  of  3  inches,  a  100  kt.  wind,  no  wave  and  a  surface 
current  of  6  kts  As  before,  the  steady  state  urogram  determined  the 
equilibrium  conditions,  which  established  the  initial  conditions  for 
the  dynamic  program  The  dynamic  program  was  then  allowed  co  run  to  a 
real  time  of  55  seconds  during  which  no  environmental  conditions  were 
changed.  The  configuration  of  the  line  cha.nged  somewhat,  but  not 
beyond  what  was  felt  were  ac—iptable  limits.  It  was  found  that  the 
angle  at  the  top  of  the  line  decreased  (which  was  anticipated  because 
of  the  change  from  equal  angle  increments  to  eauai  distance  increments 
along  the  line)  as  well  as  the  angle  at  the  bottom  of  the  line.  The 
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prograa  ran  until  the  angle  at  the  top  decreased  to  a  mimimum  value, 
and  then  was  increasing  to  when  the  con5)utations  were  stopped.  The 
line  tensions  were  re-distributed  somewhat  but  the  total  line  length 
and  the  buoy  coordinates  changed  very  little.  A  summary  of  the 
computations  is  given  in  Table  5. 

It  was  next  desired  to  subject  the  five  mooring  conditions  to 
waves  of  various  frequencies,  or  lengths.  In  each  case  the  wave 
height  to  length  ratio  was  kept  constant  at  1:15.  Thus  the  waves  were 
not  nearly  breaking  but  also  they  could  not  be  classified  as  small 
amplitude  waves.  The  dynamic  tensions  at  various  stations  along  the 
line  were  determined  and  plotted  by  the  computer.  Generally,  the 
mrximun  dynamic  tensions  occurred  at  the  buoy  or  at  the  anchor. 

Examples  of  tensions  at  the  buoy  and  at  the  anchor  are  presented  as 
Figs.  17  through  52.  Point  31A  refers  to  the  attachment  point  at  the 
buoy.  It  will  be  noticed  that  other  frequencies,  in  addition  to  the 
frequency  of  the  iropc-sed  wave,  are  present.  All  periods  of  frequencies 
in  evidence  will  be  accounted  for  later. 

The  computer  outpup  included  tension  contours  in  the  s-t  plane. 
Examples  are  shown  for  conditions  I  through  V  for  the  SCO  feet  long 
wave  as  Figs.  33  through  37.  The  value  of  each  contour  line  has  not 
been  iidicated  but  can  be  determined  by  the  reader  by  referring  to 
Figs.  17  through  32.  Tne  impoi cant  thing  in  Figs.  33  through  37  is 
to  notice  the  pattern  of  line  tensions.  The  nodes  of  zero  or  small 
tension  fluctuations  are  clearly  visible,  as  are  the  regions  of 
maximum  tension  fluctuations.  Equations  35  or  49  can  also  be  used  to 
approximate  the  locations  of  the  nodes  if  the  z  is  replaced  with  s 
That  is,  the  nodes  occur  where  cos  ws./a  =  0  .  A  comparison  of  the 


34 


results  obtained  from  the  numerical  model  with  those  from  Eqs.  33 
and  49  is  presented  in  Table  6.  Tlie  comparison  is  fairly  good  except 
at  the  highest  wave  frequencies.  Thus  the  position  of  the  nodes  o,an 
be  estimated  from  Eq.  33,  regardless  of  the  line  curvature.  This 
introduces  the  possibility  of  determining  the  best  position  for  locat¬ 
ing  sensitive  instrumentation,  at  least  with  respect  to  line  position. 
That  is,  the  positions  of  ainimum  displacement  variation  can  be 
determined  from  Eqs.  51  and  44  by  setting  sin  ^  =  0  ,  or  the 

minimum  tension  variation  can  bs  determined  by  setting  cos  —  =  0  or 

Bi 

the  minimum  line  velocities  can  be  determined  from  a  contour  map  of 
velocities  or  by  taking  the  total  time  derivatives  of  Eq.  31  or  44. 

A  conqjlete  summary  of  the  runs  made  with  waves  for  conditions 
I  through  y  is  given  in  Tables  ^  through  11.  Included  in  the  Tables 
is  the  value  of  the  damping  coefficient,  Q  ,  used,  the  range,  or 
double  amplitude,  oi  dynamic  tension,  the  mean  value  of  line  tensions 
at  the  end  of  the  mn,  the  periods  of  the  frequencies  evident  in  the 
output  records  and  a  calculated  line  tension  based  on  Eq.  33  for  which 
the  length  cf  the  line  was  taken  as  the  "stretched"  or  "in-place" 
length  at  time  ecual  to  zero,  and  the  modulus  of  elasticity  was  taken 
as  that  corresponding  to  the  tension  in  the  top  of  the  line  at  time 
equal  to  zero. 

It  is  de.sirable  to  be  able  to  display  and  compa-e  mooring  line 
tensions  in  a  way  that  takes  into  account  the  difference  in  line 
diameter  and  length,  the  differences  in  wave  height  and  frequency  end 
the  resonant  frequencies.  Thus  it  is  desirable  to  normalize  the 
tensions  and  corresponding  frequencies  to  compare  the  results  cf 
different  types  of  moorings.  The  normalizing  factor  for  tensions  for 
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this  stu'iy  vfas  based  on  thi  total  saturated  n.’. ;s  of  the  line 
(excluding  the  added  mass)  and  a  characteristic  wave  acceleration. 

Thus 

Normalizing  force  =  y  L  .H  .  (50) 

Since  eq.  33  can  predict  the  nodes  in  the  line  fairly  well  for  the 

moorings  presented  here;  it  was  felt  that  Eq.  34  may  suffice  for 

predicting,  at  least  approximately,  the  resonant  frequencies  with 

regard  to  line  tensions  despite  the  presence  of  considerable  curvature 

in  the  line.  Thus  the  first  modal  frequency  predicted  by  Eq.  34  was 

used  as  the  normalizing  frequency.  Conveniently,  the  first  resonant 

condition  should  occur  at  f/f^  =  1.0  and  the  subsequent  higher  mode 

resonant  conditions  occur  at  f/f.  =  2.0,  3.0,  4.0,  etc.  The  resonant 

frequencies  of  interest  for  the  five  conditions  are  presented  in 

Table  12.  The  resulting  normalized  frequency  response  curves  for 

conditions  I  through  V  are  presented  as  Fig.  38  for  mooring  line 

tension  at.  the  buoy  and  Fig.  39  for  mooring  line  tension  at  the  anchor. 

Admittedly,  the  data  is  sparse,  but  the  trend  for  moorings  considered 

here  is  clear.  That  is,  the  greatest  dynamic  response,  especially  at 

the  resonant  conditions,  occurs  for  the  higher  wind  conditions.  ITie 

next  most  influential  parameter  is  the  line  scope.  It  has  already 

been  seen,  and  is  common  knowledge,  that  scope  is  a  most  important 

parameter  for  steady  state  loads.  The  responses  at  higher  frequencies 

tend  to  be  damped  out  due  to  hydrodynamic  action.  So  data  was 

obtained  for  f/f  =  1.0  because  such  a  condition  would  be  due  to 
n 

waves  considerably  longer  than  1000  feet  and  such  long  waves  were 
uneconomical  to  consider.  However,  little  energy  exists  for  waves 
longer  than  1500  to  2000  feet  in  roost  wave  spectra. 
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With  ssveval  curves  like  Figs.  58  and  39  for  a  coaiplete  rang©  of 
water  depths,  scopes,  line  dia:neters  and  t>'pes,  it  may  be  possible  to 
generate  a  non-dimensional  e^3piricai  solution  for  dynamic  line  tensions 
for  any  homo<;eneoiiS  one  point  Tsooring.  Or  perhaps  analytical  transfer 
functions  lased  on  linear  systems,  which  are  being  developed  by  others, 
can  be  compared  to  results  like  Figs.  33  and  39,  which  consider  most 
of  the  non-linearities  in  the  problem.  This  subject  has  been  left 
for  future  study. 

in  order  to  test  the  program  for  quits  shallow  water  conditions 
and  to  gain  whatever  information  possible  with  just  two  runs,  condition 
VI  for  a  depth  of  150  feet  and  condition  VI i  for  a  depth  of  1000  feet 
were  investigated  for  a  wave  6Q0  feet  long  and  40  feet  high. 

For  the  depth  of  150  feet  the  sequential  positions  of  the  line 
are  shown  in  F.-g.  40  from  time  zero  to  13  .seconds.  At  about  time 
equal  to  14  seconds  the  line  tension  became  negative  and  the  program 
stopped  because  a  square  root  of  the  negative  tension  occurs.  The 
line  tension,  which  was  nearly  consta.nt  along  the  line  at  any  time  is 
shown  in  Fig,  41. 

For  the  depth  of  1000  feet  the  line  ''emained  nearly  straight  for 
all  motion.  Tne  line  tensions  again  were  nearly  constant  along  the 
line  at  any  time.  The  fluctuations  in  the  tension  at  the  buoy  and  at 
the  anchor  are  shown  in  Fig.  42.  It  should  be  recalled  for  both 
conditions  VI  and  VII  that  the  wind  velocity  was  50  kts. 


SlMdARY  AND  CONSLUSIONS 


A  review  of  the  basic  equations  for  the  numerical  modeling  of  the 
mooring  line  motion  based  on  the  method  of  characteristics  and  the 
motion  of  the  buoy  has  been  made.  Solutions  for  a  homogeneous  mooring 
line  have  been  presented  whi'-.h  include  hyateretic  damping  and  hydro- 
dynamic  damping.  In  addition,  a  summary  of  analytical  solutions  for 
a  straight,  taut  laooring  line  has  been  presented  and  it  was  attempted 
to  apply  the  procedures  developed  to  the  single  point  curved,  or  less 
taut,  nsoorings.  It  was  found  that  this  approximate  procedure  did 
poorly  for  estimating  mooring  line  tension  but  did  quite  well  in 
estimating  the  positions  of  the  nodes  of  tension  and  in  predicting  the 
natural  resonant  frequencies  of  the  line. 

A  new  way  of  presenting  frequency  response  curves  for  taut  or 
slack  mooring  lines  subjected  to  waves  was  suggested.  The  normalizing 
procedure  makes  it  possible  to  show  the  frequency  response  for  lines 
of  different  lengths,  diameters  and  materials  on  the  same  scale.  Thus 
the  influence  of  scope  and  wind  magnitude  on  dynamic  mooring  line 
tension  was  presented  while  normalizing  the  influence  of  line  diameter. 
It  was  found  that  smaller  scope  and/or  higher  wind  increased  the 
dynamic  tensions  drastically,  as  well  as  the  steady  state,  or  mean 
line  tension. 

Two  uns  were  presented  for  shallow  water.  It  was  seen,  that  a 
large  wave  in  ISO  feet  of  water  wouid  have  been  disastrous  to  a  2  inch 
difuaeter  line.  However,  the  same  wave  in  1000  feet  of  water  with  the 
same  line  and  scope  produced  reasonable  dynamic  tensions.  It  was 
seen  that  the  tension  was  nearly  constant  along  the  line  at  any  time. 
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Probably  a  static  solution  approach  would  be  adequate  wherein  the 
line  tension  would  be  dctenained  based  on  the  seouentiai  positions 
of  the  buoy  on  the  wave  and  tha  subsaquent  change  in  scope  without 
regxjrd  to  accelerations,  However,  this  procedure  would  be  very  inade¬ 
quate  for  deep  water  worings. 

A  presentation  of  the  water  sheave  problem  was  made  which 
illustrated  an  additional  application  tor  cht  numerical  program. 
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TABLE  ]  -  -  COMPAlilSON  8ETrfEE.<i  TOEOkY  AND  THE  N'JMERICAL  MODE!,  POR 
lONGI  UDINAL  DYNAMIC  TENSIONS  IN  A  '^TPhIGHT  LINE, 


Depth  ox 

L.’nc  Length 

Variation  in  Lire  Tension 
at  the  Ai>chor  (lbs) 

Variation  in 
at  the  Top  of 
(lbs) 

Line  Tension 
the  Line 

Num.  Model  Eq.  33 

Num.  f<*‘.tdei 

Eq.  .33 

5,000 

L  ,380 

1 . 330 

296 

r> 

10,000 

3,250 

3,200 

2,990 

2,940 

15,000 

1,640 

1,600 

990 

932 

20..00C 

l,7fn 

1.750 

1,1/0 

1,172 

Line  diameter  =  1.5  inches 
Modulus  of  Elasticity  -  13.1  x 

10^  psf 

Line  density  =  0.0243S  slugs/ft  fsame  as  sea  water) 
Anpiitude  of  forcing  action  =  15  ft 
Frequency  of  forcing  action  =  0.7  rad/ sec 


TABLE  2.  -  COMPARISION  BETWEEN  NATUBj\LLY  DAMPED  A.ND  FORCEFULLY  DAMPED 
LONGITUDINAL  VIBRATIONS  FOR  A  DEPTH  OF  20,000  FEET. 


Value  of  the 
Damping  Coeff . , 

Q 


Variation  in  Line  Tens  Variation  in  Line  T^^asion 

at  the  Anchor  (lbs)  at  the  top  of  the  Line 


Num.  Model  Eq.  33  .Num.  Model  Eq.  33 


0 

4x10^ 

(D 


3,237  ®  3,677 

2,085  N.A.  2,348 


N.A. 


Same  line  and  forcing  characteristics  as  for  Table  1,  e**cept. 
Frequency  =  0.6103  rad/ sec 
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TABLE  3.  -  SIM4ARY  OF  STEADY  LOAD  LONDITJONS  FOR  THE  FIVE  MOORINGS  IN 
20,000  ft.  Time  «  0 


Condition  Kind  Line  "Stretched"  Tension  at  Tension  at 

No.  Vel.  Dia.  Scope  Scope  Buoy  (lbs)  Anchor  (lbs) 


(Kts) 

(in) 

S  UIT 

%  UIT 

I 

150 

3... 

0.80 

1.05 

99,804 

.47 

86,614 

.41 

n 

150 

2.5 

0.87 

l.il 

47,777 

.44 

41,077 

.38 

III 

50 

2.5 

0.87 

1.05 

16,177 
.  15 

9.4:: 

.09 

IV 

150 

2.0 

1.18 

l.EC 

27,699 

.40 

2  3. .596 
.34 

V 

50 

2.0 

1.16 

1.54 

6,983 

.10 

2,635 

.04 

Note:  Scope  =  Length  of  line  on  land  t  water  depth 

"Stretched"  scope  =  Length  of  line  under  load,  in  place  i  water 
depth 
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TABLF  4.  -  SIM4ARY  OF  STEADY  LOAD  CONDITIONS  FOR  TflE  TWO  MOORINGS  IN 
SHALLOW  WATER.  Time  =  0  Wind  Velocity  =  50  kts 


Condition 

No, 

Water 

Depth 

(r-t) 

Line 

Dia. 

(ins) 

tt 

Scope 

Stretched" 

Scope 

Tension  at 
Buoy  (lbs) 

VUH 

Tension  at 
Anchor  (lbs) 
CUH 

VI 

ISO 

2 

1  32 

1,50 

2,021 

2,001 

.03 

.03 

VII 

1,000 

2 

1.32 

1.50 

2,564 

2,394 

.04 

.03 

TABLE  5. 

-  COMPARISON  OF  EQUILiBRIiM  CONDITIONS  BETWEEN  THE  STEADY  STATE 
PROGRAM  AND  THE  DYNAMIC  PROGRAM 

Program 

Time 

(secs 

Angle 
)  at 
Buoy 
(Rad) 

Angle 

at 

Anchor 

(Rad) 

Tension 

at 

Buoy 

(lbs) 

Tension 

at 

Anchor 

(lbs) 

Total  X 

Line 
Length 
(ft) 

-Coordinate 
of  Buoy 
(ft) 

Steady 

0 

1.089 

0.097 

19,977 

15,383 

24,527 

21,531 

State 

Dynamic 

65 

0.9C2 

0.021 

21,817 

12,141 

£4,252 

21,320 

Steady 

1.21 

4.66 

0.915 

1.265 

1.002 

1.010 

Dynamic 

TABLE  6.  - 

COMPARISON  Of 
EQ.  33 

TENSION  NODES  BETWEEN 

THE  NUMERICAL  MODEL  AND 

Cond i t ’ on 

Wave 

Length 

(ft) 

Node  Position  from 
NTimerical  Model 
(ft) 

Node  Position  from 

Eq.  33 

(ft) 

I 

500 

4.800 

6,920 

18,200 

20,600  (at  buoy) 

II 

145 

4,440 

3,680 

11,600 

11,000 

18,900 

18,400 

325 

5,900 

5,500 

17,000 

26,5Cr 

SUO 

6,240 

6,800 

19,500 

20,400 

1,000 

10,600 

9,650 

in 

145 

4,800  ? 

7,400 

10,300  ? 

12,300 

14,400  ? 

17,200 

325 

3-910 

3,690 

11,700 

11,100 

500 

5,150 

5,150 

15.100 

15,400 

630 

5,500 

4,600 

16,000 

13,800 

1,000 

7,550 

6,450 

8.700 


J5,300 


■  .r* 


HO 


TABLE  6.  -  Continued 


Wave 

Condition  Length 
(ft) 


IV  ISO 


27S 


500 


600 


V  150 


275 


460 


500 


Node  Position  from  Node  Position  from 

Nuaerical  model  Hq.  33 


(ft) 

(ft) 

4,020 

3,730 

11,200 

17,250 

18,700 

7,030 

5,050 

17,400 

15,200 

8,003 

6,800 

21,600 

2J,100 

8,350 

7,450 

23,300 

22,300 

2,510 

9,800 

7,550 

12,500 

3,400 

3,700 

10,200 

4,450  ? 

4,360 

13,100 

4,450 

4,580 

13,700 

4,450 

5,000 

10,500  ? 

15,000 

600 


48 


o 

u 


u  c 
-  o 

6  -H  4K 

c  c  c* 
X  a  « 
O  H  u 


a 

X 

o 

u 


55 

O 

I 

tj 

X 

ui 


> 

< 

3 


O 

X 


c/) 

a 

o 


o 

o 

o 


o 

o 

o 


o  o 
o  o 
00  rf 


o  o 
o  o 

O  iA 
sO 


«o  o 
o  o 


o  o 
o  o 

O  t/> 
o  ««t 


o  o 
o  o 
00  00 
•»  to 


O  00 

o  > 
(M  a% 


o  o 
O  </) 

l/> 

l/J  <N 


O  o 

r*  CO 


o  o 
o  o 

O  lO 


o  o 
o  o 
o  o 

^  vO 


o  o 
o  o 
o  tn 

00  'O 


o  o 
o  o 
o  tn 

IS.  nO 


o  o 
o  o 
r>>  00 
to  o 


s  <=> 
2  "=> 
2 

00  00 


o  o 
o  o 

00  n 

l/>  t/> 

•t 


o  o 
o  o 
oo  ^ 
o  crv 


o  o 
o  o 
o  o* 
<0  a» 
—  ^ 


o  o 

O  U) 

^  o 

»<>  hs 


o  o 
o  o 

nO  00 

LO  o 

F.4  lA 


o 

X 


o 

c 

o 


TABl.l.  9.  -  link  tensions  AND  PERIODS  FROM  WIND  AND  WAV  EXCITATION  FOR  CONDI  FION  III 
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RESONANT  FREQUENCIES  FOR  THE  VARIOUS  CONDITIONS,  BASED  ON  EQ.  34 
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FIGURES 


?IG.  1  DEFINITION  SKETOI  OF  THE  COOfiDINATE  SYSTSfii. 


Stress  (  psf  x  10"®  ) 


FIG.  3  STRESS-STRAIN  DIAGRAM- 


y//////////// 


IDEALIZED  VERTICAL 


FIG.  6  DYNAMIC  LINE  TENSION  AT  THE  BUOY.  WATER  DEPTH  »  1500  FEET. 


Morizontai  Displacement  (ft) 


FIG.  7  HORIZO.NTAL  LINE  DISPLACEMENT  FOR  A 
SINUSOIDALLY  DISTRIBUTED  INITIAL  VELOCITY. 
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FIG.  14  SEQUENTIAL  POSITIONS  OF  THE  1*3  INCH  LINE  FOR  THE  WATER  SHEAVE  PROBLEM. 


Dis1arc0  Botween  Line  Stations  Is  866. 
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TENSION  CONTOURS  IN  THE  s-t  PLANE  FOR  THE  WATER  SHEAVE  PROBLEM. 
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TK.NSION  AT  THE  '‘NCHOR.  CONDITION  I,  WAVE  LENGTH  --  500  ET.  WAVb  HEIGHT 
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(•IG.  27>  TENSION  AT  THE  BUOY.  CONDITION  II.  WAVE  LENGTH  =  500  FI.  WAVE  HEIGHT 
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FIG.  25  TENSION  AT  THE  BUOY.  CONDITION  II.  WAVE  LENGTH  “  1000 


ANCHOR.  CONDITION  11.  WAVE  LENGTH  =  1000  IT.  WAVE  HEIGHT  =  06. 
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29  TENSION  AT  THE  BUOY.  CONDITION  iV.  WAVE  LENCTM 
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Fir,.  SQ  TENSION  AT  THE  ANCHOR.  CO.NDITION  IV.  WAVE  LE.NGTH  =  500  FT.  WAVE  HEICHT 


riG.  3J  TENSION  AT  THE  BUOY.  CONDITION  V.  WAVE 


TENSION  AT  THE  ANCHOR.  CONDITION  V.  WAVE  LENC.TH  SOO  FI.  WAVE  HEIGHT 


CONDITION  II.  WAVE  LENGTH  --  500  FT.  WAVE  HEIGHT  =  33. 
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FIG.  38  ESTIMATED  FREQUENCY  RESPONSE  CURVES  FOR  LINE  TENSION 
AT  THE  BUOY.  DEPTH  =  20.000  FT. 


Condition 


ESTIMATED  FREQUENCY  RESPONSE  CURVES  FOR  LINE  TENSION 
LHOR.  DEPTH  =  20,000  FT. 
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FIG.  40  SEQUENTIAL  POSITIONS  OF  THE  MOORING  LINE.  DEPTH  = 

ISO  FT.  WAVE  LENGTH  =  600  FT.  WAVE  HEIGHT  =  40  FT.  WIND  =  50  KTS. 


Line  Tension  x  )0~^  (lbs) 


Time  (seconds) 


FIG.  41  ^fOORI^'GS  LINE  TENSION.  DEPni  =  ISO  FT.  WAVE 
LENGTH  ^  600  FT.  WAVE  HEIGHT  =  40  .n-.  WI.S'D  =  50  KTS. 


